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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEMORANDUM 1303 

RESISTANCE OF CASCADE OF AIRFOILS IN GAS 
STREAM AT SUBSONIC VELOCITY* 

By L. G. Loitsianskii 


A method of computing the resistance of a cascade of airfoils 
in a viscous compressible gas flow is described. 

The case of an incompressible gas is considered in reference 1 
and appears herein only as a simple particular case of the general 
theory of resistance of a cascade in a compressible gas. 

The investigation was restricted to subsonic velocities (that 
is, when the local velocity of sound is nowhere reached on the air- 
foil surface) because the required assumption of isentropic flow, 
that is, the absence of shock waves in any region of the motion, 
is valid only under these conditions. 

The second reason for the restriction to relatively small 
values of Mach number is the possibility under this assumption of 
applying a lift formula analogous to the well-known Joukowsky 
formula (reference 2) and of thus assigning a definite meaning to 
the term "cascade resistance" or, more accurately, the "resistance 
of an airfoil in cascade." 

The resistance formula can be derived for an isolated airfoil, 
as is known, by applying the moment-urn theorem between two parallel 
cross sections of the flow at an infinite distance upstream and down- 
stream of the airfoil. In the problem of cascade resistance, dif- 
ficulty is encountered, namely, the absence of an external potential 
flow downstream of the cascade where the boundary layers (wakes) 
from the individual airfoils merge. This essential difficulty, 
which is expressed quantitatively in the impossibility of employing 
the boundary- layer (wake) equation up to a plane at an infinite 
distance, can be circumvented by introducing the plane of merging 
of the boundary layers (wakes) and by establishing relations 
between the gas dynamic elements in this plane and in the plane 
at infinity downstream of the cascade. 


* Soprotivlenie reshetki profilei v gazovom potoke s dokriti- 
eheskimi skorostiami, Prikladnaia Matematika i Mekhanika, vol. XIII, 
no. 2, 1949 
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An essential assumption of the present investigation is that 
a small degree of nonhomogeneity of the flow exists in the section 
of the aerodynamic wake of the cascade where the boundary layers 
from the individual airfoils, considered as layers of finite thick- 
ness, merger the larger powers of the small velocity differences 
may then be neglected. The same assumption was made in the investi- 
gation of cascade resistance in an incompressible gas (reference l) 
and was subsequently confirmed experimentally. The plane of merging 
of the boundary layers is then assumed to be the control surface 
required for the application of the momentum theorem and in the case 
of the isolated airfoil is taken to be the plane at an infinite 
distance downstream of the airfoil. It is evident that when the 
relative pitch of the cascade is increased, this plane will be farther 
and farther away from the axial plane of the cascade and in the limit, 
for a relative pitch equal to infinity, that is, in the case of an 
isolated airfoil, will go to infinity. This assumption may evi- 
dently be made for cascades with moderate solidities, a case that 
corresponds in practice to turbine and compressor cascades. 

Any method of calculating the boundary layer in a compressible 
gas may be used to compute the characteristic thicknesses of the 
layer and to estimate the effect of the compressibility of the gas 
on the external flow. The solution of the proposed problem reduces 
to a straightforward and direct form that is independent of the 
method of computation. 

1. Resistance of airfoil in cascade. Joukowsky force as com- 
ponent of total force exerted by incompressible fluid on airfoil. - 
For two-dimensional flow of a real fluid, the resistance (or drag) 
of an isolated cylindrical wing of infinite span referred to unit 
length of the wing is the component of the total force exerted by 
the fluid on the wing in the direction of the velocity of the approach- 
ing flow at infinity, or, in other words, of the velocity component 
of motion of the wing in an incompressible medium. 

This definition is invalid in the case of an airfoil in a 
two-dimensional cascade, because in this more general case there is 
no unique velocity direction at infinity upstream and downstream of 
the cascade and there are no considerations by which preference is 
to be given to any particular direction for determining in this 
direction the resistance component of the total force acting on the 
wing. In this case, the problem is to determine what may be termed 
resistance. 


An Isolated ■'/ring of finite span is now considered. In this case, 
as also In the case of an airfoil in cascade, for each section of the 
wing, in view of the presence of vortex systems (films) shed from the 
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wing and passing downstream to infinity, two velocities different in 
magnitude and direction exist at infinity upstream and downstream 
of the wing. For ideal flow about a wing of finite span in accordance 
with the scheme of lifting lines, the total pressure force of the 
flow at a given section of the wing is known to be perpend! cular to 
the velocity of the flow at, the corresponding point of the section 
under consideration on the lifting line. This velocity, which repre- 
sents half the vector sum of the velocities at infinity upstream and 
downstream of the wing, is assumed at the section considered as the 
effective velocity of flow; the angle between the chord of the wing 
section and the direction of the effective velocity is considered 
as the effective angle of attack, and so forth. 

For a two-dimensional infinite cascade of airfoils, a similar 
assumption is made with the difference only that in the theory of 
the wing of finite span the effective velocity differs slightly from 
the velocity of the approaching flow; whereas in the case of the 
cascade the jump is of the same order as the geometrical angle 
of attack. 


In the aerodynamics of a wing of finite span, the profile drag 
is the difference between the head resistance, which is represented 
by the component of the total force exerted by the real (viscous) 
gas flow in the direction of the velocity at infinity upstream of 
the wing, and the induced drag, which is the component in the same 
direction of the effective lift force. 

For a small difference between the directions of the effective 
velocity and the velocity of the approaching flow, this definition 
of the profile drag of a wing section differs by small terms of 
higher order from the true profile drag, strictly defined as the 
vector difference between the total force exerted by the real flow 
on the wing section and the effective lift force for a real fluid. 

In the case of the two-dimensional cascade, it is natural to 
assume for the profile drag R' the difference between the vector 
of the total force R (fig. l) and the Joukowsky force R- (in 

«J 

the terminology of reference l) which for an incompressible gas 
is given by 


R. = 
J 


pv r 

m 


( 1 . 1 ) 
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acting in the direction perpendicular to the fictitious velocity 
at infinity V m determined as 


* 

where (u.^, v^) and (u^, v^) are the vector velocities 

at infinity upstream and downstream of the cascade, p is the den- 
sity of the fluid, and T is the circulation determined by the 
equation 


r “ ( v 2oo “ C 1 * 3 ) 

where t is the pitch of the cascade. 

Introducing the concept of the vector pitch t, which is equal 
in length to the magnitude of the pitch t and directed at right 
angles to the axis of the cascade downstream of the flow, gives 
the Joukowsky force by the following vector equation (reference 2): 


= pV m x (t x 


m 


V 


(1.4) 


where the following vector 


V, 


V - V 
* 2oo v l® 


(1.5) 


gives the vector change of velocity produced by the cascade. The 
preceding formulas are valid not only for the flow of an ideal incom- 
pressible fluid but also for a viscous fluid. 

The profile drag R* is as follows (reference l): 


R* = p’t (1.6) 

where p' is the pressure loss in the cascade determined by the 
equation 
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P' = (Pl» + | C'^J) - (P2« - 1 P V 2») 

- ( P l. + I P \J) - ( P 2. + I pV 2» 2 ) 

The total force R is equal to the sum 


(1-7) 


R = Rj + R’ = pV m X (t X V d ) + p't (1.8) 

2. Resistance of two-dimensional cascade in real gas flow at 
subsonic velocities . - The expression for the total force R of 
the interaction of the flow with a two-dimensional cascade at sub- 
sonic velocities may be represented in the following form (reference 2) 


B - - PfJ* + < v i„ • - p 2» < v a» • f 2 - 1 ) 

where p, , p n and p 0 , p 0 are the pressures and densities 

JSX> L.CD oOO 

upstream and downstream of the cascade, respectively. The equivalent 
expressions 


p loo V l«* ’ t = p 2e* V 2o> • t 


(2.2) 


evidently express the rate of mass flow per second through the 
section of the flow parallel to the axis of the cascade and equal 
in length to the pitch. 

As was shown (reference 2) also in the case of a compressible 
gas for Mach numbers not too near unity, the lift force of an airfoil 
in cascade in an ideal gas flow may be represented in the form of 
equation (1.4), provided that for the density p is taken the arith- 
metical mean density equal to 

■°m = 2 ( p lc*> + p 2oJ ( 2 *3) 


The following approximate expression of the Bernoulli theorem 
is employed: 


p. - p Q = p V • V = i p (V 2 - V, d ) 
loo 2 od m d 2 m 2oo loo 


(2.4) 
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This equation is valid with an accuracy to tenth parts of the square 
of the difference of squares of the Mach numbers at infinity upstream 
and downstream of the cascade. 

In the case of the real (compressible and viscous) gas, 

•p, - p, = - p (V, 2 - V, 2 ) + p' (2.5 

~ loo 2ot> 2 m 2oo loo r v 

where p f characterizes the losses in the cascade due to the internal 
friction in the gas; an equation may be obtained (reference 2) analo- 
gous to equation (1.8) 


E=R. + E' = pV X (t X Vj + p'- 
j mm d 


( 2 . 6 ) 


where p' is determined by the expression 


p ' - P l» p 2oo " \ °m (V 2. 2 - V 1» Z) 


(2.7) 


The problem of determining the profile drag force R' equal to 


R' = p't (2.8) 

thus reduces to finding the losses p' which depend on the shape 
of the airfoil in the cascade and the character of the flow about 
the airfoil. 

3 . Introduction of intermediate plane ; relation between gas 
dynamic elements in this plane and corresponding values at infinite 

distance from cascade . - In addition to the planes lec and 2 os 
that were employed in the analysis of the incompressible fluid, 
(reference 1) an intermediate plane 2 is introduced for the compress- 
ible gas (fig. 2); plane 2 is located where the boundary layers 
(wakes) from the individual airfoils merge. The hydrodynamic and 
thermal boundary layers in the wake downstream of the cascade are 
hereinafter assumed to have the same pattern. 

The following assumptions with regard to the motion of the gas 
near plane 2 are necessary: By definition of the position of plane 2, 

no individual boundary layers exist in the flow downstream of this 
plane; the aerodynamic and thermal wakes of the airfoils are, however. 
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maintained and depressions in the velocity or total-pressure curves 
and also depressions or peaks in the temperature curves result. 

A fundamental property of the boundary layer is that the pressure 
transverse to the wake is the same at all points of a given normal 
section of the wake; that is, no pressure drop in the distribution 
curve occurs in this section. The pressure along the wake changes 
sharply in the immediate neighborhood of the trailing edge of the 
airfoil and is gradually equalized as the distance from the trailing 
edge is increased. 

Two sections of the wake are passed through the point of inter- 
section of plane 2 with the axis of the wake; one section lies in 
plane 2 and includes the y-axis (fig. 3), and the second section 
lies in plane 2' normal to the axis of the wake and includes the 
y ’-axis. 


The following magnitudes are introduced 

u 2 “ u 


A 

u 


-if 


y 0 +t 


u.- 


.y 


0 


2 


dy 


-if 


y o +t 


P2 - P 


dy 


-if 


y 0 +t 


v 2 " 


v 


v t 


A = 


f 


'0 

y 0 +t 


v 2 


dy 


y 0 


P2 - P 
P2 


dy 


t 


f 


Yr 


y 0 +t 


v 2 - V 


dy 


y o +t 


^T 


-if 


To - T 


T 2 


dy 


’0 


r (3-i) 


j 
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which characterize the mean relative deviations of the hydrodynamic 
elements of the flow at the points of section 2 of the wake from 
the values of these elements at the boundaries of the wake at the 
points of intersection of the boundary layers. 

Section 2 will be assumed at such distance from the cascade 
that the differences u 2 - u, . . . , and also their mean relative 

values A u > ... may be considered small magnitudes, the higher 

powers and the products of which may be neglected. Moreover, the 
velocities at different points of section 2 are assumed parallel 
and in a general direction coinciding with that of the velocity at 
infinity behind the cascade. It follows at once that 


A u= A ,= A V 


(3.2) 


Comparison with analogous mean relative deviations in section 2* 
gives the magnitudes 


A u’ 


4J 1 

^0 


y '+-t 


u ’ - u* 


u.- 


dy 1 , . . . (t* = t cos P 2oo ) 


(3.3) 


In the subsequent discussion, it will be assumed that, for a 
sufficient distance of planes 2 and 2’ from the axial plane of the 
cascade, all the magnitudes (3.3) and so forth are correspondingly 
equal to the magnitudes (3.1); that is. 


A u * -A u Ap' -Ap • • • (3.4) 

This additional assumption may be justified as a consequence 
of the assumption of a small degree of variation of the gas dynamic 
elements near plane 2 and behind it downstream of the flow. 

In accordance with the fundamental property of the wake /\^ ’ =0, 
the following equation may be obtained: 


Ap = 0 


(3.5) 



NACA TM 1303 


9 


Because of the smallness of the magnitudes A a , A p , . . . , the 

gas dynamic magnitudes in the intermediate plane 2 are easily shown to 
he connected with the corresponding values of these magnitudes in the 
plane 2oo by relations that are analogous to the case of the incom- 
pressible gas. 


For this purpose, a segment of a flow tube is assumed between sec- 
tions 2 and 2«*>, where a length equal to the pitch t is taken for the 
transverse dimension of the tube in the direction parallel to the axis 
of the cascade. 


Application of the theorem of the conservation of mass then yields 

y 0 +t y 0 +t 

J pu dy — j [p 2 - (p 2 - p)] [u 2 - (u 2 - u)]dy = P^u^t 
y ° y 0 

Expanding the brackets and neglecting the product (p 2 - p) (u 2 - u) 
as a small quantity of higher order gives the following equation: 


y 0 +t 


[p 2 u 2 - u 2 (p 2 - p) - p 2 (u 2 - u)]dy = P^u^l 


From this expression, the following relation is obtained in the 
notation of equation (3.l): 


P 2 U 2( 1 " A p “ ~ p 2oo U 2oo 

or, with the same degree of accuracy, 

P 2 U 2 = P2oo U 2oo^ 1 + A p + V ( 3 * 6 ) 

The momentum theorem in the projection on the x-axis applied to 
the same segment of the flow tube gives 

r 0+t r 0+t 2 

J pdy +J pu2dy - p 2o t - Pg^ug^t = 0 

yo y 0 
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This equation may be written in the form 
y 0 +t y 0 +t 

J &2 - C p 2 - P IK + J [ p 2 - ( p 2 - p 0 C U 2 - C u 2 - u H dy 

y 0 y 0 


p 2 c»? + P 2 od U 2 oo t 


If the smallness of the differences P 2 - P> P 2 “ P> and - u 2 
is taken into account, the following expression may be obtained: 


- u 


p 2^ “ V + P 2 U 2 2 ^ ' A p - ^u) “ P 2oo + p 2» U 2o 0 2 ( 3 * ? ) 

With the aid of equations (3.5) and (3.6) and the same approxima- 
tion, the following equation may be written: 


P2 + P2oo u 2 0 o u 2( 1 - An) = P2oo + P2oo u 2os^ 


(3.8) 


The mementum theorem is now applied in the projection to the 
y-axis, which gives 


X 


y o +t 


puv dy == P 2 o u 2 o v 2oo t 


y 0 


or 

1 

t 


y 0 +t 


J [P 2 " (P2 - P>] [ u 2 - < u 2 “ u fl " ( v 2 ~ v 3 dy = P2a u 2x V 2o 0 


y 0 


Rejection of small terms of higher order leads to the equation 


PZ^^C 1 - A d " = P2od u 2oo v 2c 
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or, according to equation (3.6), 


v 2 (l ~ \) = 




V 2 " 


- + A v> J 


(3.9) 


The assumption of parallel directions of the velocity vectors in 
sections 2 and 2®, with the aid of equation (3.2), yields 


u 2 - u 2»h + V 

" VC 1 + V - V 2 „(l + A u ) 
Equation (3.8) then gives 


(3.10) 


P 2 P 2od 


( 3 . 11 ) 


On the basis of equation (3.10), there also follows from equa- 
tion (3.6) 


P 2 = P 2 J 1 + A p ) 

Finally, from the Clapeyron equation, 
n Pp - (Pp - P) 


( 3 . 12 ) 


p P 2 - (P 2 " P) 


= RT = R[T 2 - (T g - T)] 


or, when the smallness of the differences is accounted for. 


P2 

P2 




P 2 ~ P P2 ~ P 




P2 p 2 

From this equation, - A^ = is obtained by integration, or. 


A p 


(3.13) 
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Conversely, the same Clapeyron equation In planes 2 and 2w 
yields, "by equation (3.1l) 

R(T 0 - T ) = i-2£P - ^ -Zg. £2® ^ 

2°° 2 p 2< » P2 p 2oc p 2 p 2o= p 

or, by equation (3.13), 


t 2 » - T a - T a» A p ’ - 

that is, 

T 2 = T 2 d 1 + V (3 * 14) 

4. Relation of fictitious wake thicknesses to magnitudes A and 
A'. Expression of profile drag in terms of fictitious wake thick - 
nesses . - The momentum equation in the wake behind the airfoil of 
the cascade will now be employed; the equation contains the following 
fictitious wake thicknesses defined (reference 3) as integrals over 
section 2': displacement thickness 6 2 * and loss -of -momentum thick- 

ness 6 2 **, 



When these thicknesses are connected with the magnitudes A, 
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-r (■ • 


C p 2 ~ (P2 ~ pfl C v 2 ~ ( v 2 ~ v 0 


vf 


y 0 ,+t 


yo’ 


P2 - p 
P2 


P 2 V 2 


■} 




dy 1 


dy’ + 


f 


y 0 ,+t ' v 2 - V 


V 2 


dy 1 


> (4.2) 


= t'(A ' + %') = (A D + Ajt COS P 


2oo J 


* * 


■X 


yo'+t 


° [p 2 - (P2 - p)] [ V 2 - ( V 2 - V)] V 2 - V >. 


v 


P 2 V 2 


dy 1 


f (4.3) 


t'Ay* - tAjj COS p 2oo 




The profile drag will now he determined; the magnitude p' must 
first he found. In equation (2.7), p* is expressed as a small dif- 
ference between two large magnitudes and is therefore unsuitable 
either for experimental or for approximate theoretical determination 
of p’. In order to eliminate this defect, equation (2.7) is rewrit- 
ten in the form 


p ' “ p loo " p 2 t» ~ 4 (Pioo + p 2 co) ( V 2<» 2 " V loo 2 ) ( 4 *4) 

and the flow is considered between section lcc and the limits of the 
boundary layers that merge in plane 2. In this entire region, the 
flow is nonvortical so that the Bernoulli theorem may be applied 
without the additional term that accounts for such losses. The fol- 
lowing expression similar to equation (2.4) may then be written: 


p lco “ p 2 “ I ( p lt» + p 2^ ( V 2 2 " V lo= 2 ) 


(4.5) 
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Because by equation (3.1l) p 2 = p 2oo , the following expression 
is obtained when equation (4.5) is compared with equation (4.4) : 

*’ ‘ t 2 > - i (p„ + oa.) <V, 2 - v *) 


When V 2 and p 2 are replaced by their expressions in terms of 
V 2 o,, and Pgoo^ then according to equations (3.10) and (3.12) , 

p' - I (p to + p a» + p axA p ! [v e „ 2 (i + ! - vj 2 ] - 


4 + p aJ < v 2„ 2 - v i„ 2 ) 


- | Cpl» + p aJ v aA + j p a.( v 2» 2 - O A f 

or, by equations (4.2) and (4.3), 

1 


V* 


p' = p m v ? 2 + x Po (Vp 2 - V-, 2 ) — 

E t cos (3^ 4 ^°° 2 » t 


So* - 60** 


cos 


% 


(4.6) 


The following magnitude is now introduced 1 


60* A.. + A n 

Ho = — = £ = 1 + __P 


ft ** 
°2 


Ai 


% 


(4.7) 


which is for the case of the motion of an incompressible gas; the 
fo ll owing simple equation is then obtained: 


P' - [ p m ^ 2 + 1 (% - 1) ezJvzJ - v lre 2 )] 


6 2** 


t cos P 2oa 


The formula for profile drag is immediately obtained from 
equation (4.8), 


(4.8) 


R’ = P't = [p^ 2 + i (H 2 - 1) P - vj*)] _!a£- (4.9) 
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From this expression, the profile-drag formulas for a cascade in 
an incompressible viscous fluid are obtained as particular cases and 
for the isolated airfoil in the general case. 

For the case of a cascade in an incompressible fluid 
(p =* constant) , ■ 



and equations (4.8) or (4.9) are converted into 

p, _ ** 

t cos 0^ 

R . = pvacfoao, ** 

cos 02oo 

which are identical to equations (2.12) of reference 1. 

For an isolated airfoil in a viscous compressible fluid. 


-\ 


J 


(4.10) 


P lo> P 2oo ~ p m " P °° 


V. = V„ = V 

loo 2oo oo 


* 20 ," 0 


Moreover, plane 2 extends to infinity, so that 


R' = p v ^8 ** 

^ Nx V OD u oo 


(4.11) 



16 


NACA I'M 1303 


Equation (4.1l) is the well-known formula of the resistance theory for 
an isolated airfoil. 

The losses and the profile drag of the cascade are expressed by 
equations (4.8) and (4.9) in terms of known elements at infinity 
ahead of and behind the cascade and in terms of the elements Hg and 

5g** referred to plane 2, the position of which remains unknown, 

because up to the present no reliable theory of the turbulent wake 
exists. 

A formula will now be obtained for the profile resistance of the 
cascade; by the theory of the boundary layer at the airfoil, this 
expression makes possible the computation of the resistance of the 
cascade, and the dependence of the magnitudes Hg and &g* * just 

mentioned on the elements of the boundary layer at the rear edge of 
the airfoil of the cascade can therefore be determined. 

5. Establishment of relations between wake elements in sec- 
tion 2 and boundary- layer elements in sections at trailing edge. _ 

A generalization is given herein of the known device of setting up 
relations between the elements of the boundary layer at the trailing 
edge of the airfoil and in the wake behind it at infinity, as proposed 
for the case of the isolated airfoil in the incompressible fluid by 
reference 4. 

In this generalization, for the case of the cascade the section 
at the trailing edge is connected not with the plane at an infinite 
distance downstream of the flow but with plane 2 of the merging of 
the boundary layers or, more accurately, with plane 2 1 inclined to 
it by the angle Pgw Moreover, the generalization requires passing 
to the compressible case. 

The momentum equation for the wake behind a body may easily be 
derived from the general equations of the plane boundary layer in a 
compressible gas 


pv P 


5v„ 5V„ 

— + P V n — 


5s 


5n 


dp 5 t 
ds 5n 


a(pv s ) a( P v n ) 

5s + 5n ~ 


-\ 

> 

J 


(5.1) 
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where for the longitudinal (coordiante s) and transverse (coordinate n) 
projections of the velocities, the symbols V s and V n are used in 

contrast to the velocity projections u and v connected with the 
axes Ox and Oy; p is the local density, t the friction stress, and 
p the pressure on the outer boundary of the layer. 

By rewriting the system (5.1) in the following form, according 
to the second of equations (5.l) and the general Bernoulli equation. 


d 

Ss 


(pv s v s ) 



= pv s 


+ ar 

ds c)n 


5(pV s ) S( P v n ) 

Ss + Bn ~ 



(5.2) 


where p and V s denote the density and the longitudinal velocity 
at the outer limit of the boundary layer. Both sides of the second 
equation are then multiplied by V s to yield 


d , _ , 3 , - , dV_ 

S (pT - T ») + £ - pV = if “ 0 

The first of equations (5.2) is then subtracted term by term 
from the equation just obtained; the resulting equation is then inte- 
grated, which gives 




+ (pv s - P v s ) — S. 

ds 


5t 

dn 


along the normal to the section of the wake, which is considered either 
infinite in the usual sense of the theory of asymptotic boundary layer 
or finite, as is assumed in the. theory of the finite thickness layer. 

In either case, the following relation holds: 


+°°,8 — + 03 , S 

if pv.c*. - v s> ^ + Sr r c^ s - ^ - ° 

-oo, 8 -o,, 5 

The following expression is then obtained: 
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By expanding the parentheses and introducing the notation of 
reference 4, 


5* 



iBoj 5 



P^S 

pv s 




(5.3) 


the required momentum equation is finally obtained. 


ds \v ds 


p ds/ 


.* * 


1 dV 
V ds 


8 * 


= 0 


(5.4) 


In such form, the momentum equation for the compressible gas 
differs from the corresponding equation for the incompressible gas 
only in the term "p dp/ds (and, of course, in the definitions of 
the magnitudes 8* and 5**). If the momentum equation for the 
incompressible gas is considered for the case of axial symmetrical 

motion, the term p -1 dp/ds, which expresses the effect of the variable 
density of the gas, may be taken equivalent to the term that takes 
into account the transverse curvature. 


In addition to the momentum equation, the heat equation is con- 
sidered; it can be easily set up by a method analogous to the preceding 
method from the known heat equation of the boundary layer. 


( 


pV s ^ + pV n ^ 


:)(**«■ 


dq 

3E 


where a = pCp/k is the Prandtl number. 


(5.5) 
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The value of q is given in the case of the laminar boundary 
layer; equation (5.5) holds also for the turbulent layer, but in this 
case q would be expressed in a different form. 

The so-called temperature of adiabatic stagnation T * is now 
considered; it is given by 


T* 



(5.6) 


By means of the continuity equation, the following system of 
equations may be set up: 


(pV s T*) + S (pV n T*) = -JL- M 
ds on Jc-r. on 




^ (PV S T*) + ^ ( P V U T*) = 0 


J 


(5.7) 


In the second equation of the system, the stagnation temperature 
T* ' at the outer limit of the boundary layer, which is constant 
(because the external flow is isentropic), is taken under the sign of 
the derivatives in the continuity equation. 


Subtracting one equation of the system (5.7) from the other and 
successively integrating over the cross section of the wake gives 


+Oo, 5 

PV S (T*-T*)dn:=0 

o, & 




+ 00,5 

J"' pV s (T* - T*)dn = constant 


(5.8) 


“ 03 , <■ 




The fictitious thickness of the wake is now introduced 
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which may he termed the thickness of the energy loss; equation (5.8) 
may then he rewritten as 


pVT*0 = constant 


(5.10) 


Equation (5.4) is again considered. After each side is divided 
hy &**, the expression is integrated along the wake from section k 
at the trailing edge of the airfoil to plane 2' , previously defined. 
The result is 


In 


/2s!!) =ln /£>!£). ^ 

\6 k *V Wz 2 / « k) 


5* d In Vc 


;** 


ds 


ds 


(5.11) 


The notation of equation (4.7) is used for the ratio of the ficti- 
tious wake thicknesses , 



(5.12) 


and it is noted that equation (5.1l) is integrated to completion 
if the magnitude H is replaced hy some average value; for example, 
the following relation may he set up: 


E “ H op = I < H 2 + s k) 


(5.13) 


By this simplification, the following expression is immediately 
obtained: 


In 




i (Ho + H k ) In -2 


or finally, 



Pk 

P2 



+ | ( H 2 +H k) 


(5.14) 


This equation connects an< 3- hut does not explicitly 

contain the exponent on the right contains the magnitude Hg, 

which is equal to the ratio From equations (4.8) and (4.9) 

previously derived, equation (5.14) serves as one of the equations 
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for expressing the two unknowns 5 g ** and Hg entering In the 
equations for the losses in the cascade and the resistance in terms of 
the elements of the houndary layer at the trailing edge of the airfoil. 

The second equation is obtained by use of equation (5.10), which 
may be rewritten as follows : 

P2V 2 T 2 * 0 2 = Pk V k T kX 

or, because of the isentropic character of the motion outside the 
boundary layer, Tg* = the expression then becomes 


P 2 v 2 e 2 = Pk v k 0 k (5.15) 

In this equation a new unknown quantity 0g appears to enter; 

because of the small degree of nonhomogeneity of the fields of hydro- 
dynamic elements in planes 2 or 2 1 , however, this term can actually 
be expressed in terms of the previous unknowns. When the small degree 
of nonhomogeneity and the formulas relating the elements in planes 2 
and 2on (derived in section 3) are accounted for, equation (5.14) and 
then equation (5.15) are transformed. By equation (5.14), 

S 2 ** p k r V k “I 2 '4( H 2 +H k) ^ 

57* * f ^ 1 + v Ua 1 + vJ 



Because in this section everything will be expressed in terms of 
the unknowns Sg** and Hg, equations (4.2) and (4.3) are applied; 

the following expression is then obtained: 
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*£ = r . . u 2 £! . 

5 k** P2oc \VaJ L 2 f 




or 


a ** _ p, ** p k / v k \ 

2 k P2 oA V 2cJ 


^[ 1 -( 1 + | H k + l H 2 )£| l *] 


r ( 5 . 17 ) 


For a first approximation, the subtrahend in the brackets on the 
right may be neglected in comparison to unity to obtain 


5 **- 6 **0^/1^ 

k P2 oAW 


2 ^k4 H 2 , 


(5.18) 


The second fundamental equation (5.15) is similarly transformed. 
With the chosen degree of accuracy. 


+° 9> 6 „ . 

T* \ 

+ 03,5 

r p s ti - 


“I r- P 2 V 2 \ 

V/ 



- 03,5 

+ODJ 

- T* ) 

+ 03,5 

f 02 V 2< T 2 

: 

J R P 2 V 2 T 

* 

2 

- 03,0 


- 03,6 


P V s (T 2 * - T*) 

P' 2 V 2 T 2 * 


VLl^dn 

V 


dn 
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Therefore, 


+oo,5 


+oq.5 


f J5z_ 

J c To + 


V? , 2 /2JCn - Vsf/gJCjB 


fog - T)dn ^ P 

,8 "2 + V 2 2/2Jc p J^ b T 2 + V 2 2/2J Ci 


dn 


+to,5 


= — 1 P dn + 

1 + V 2 2/2Je p I 2 *2 


+oo, 6 




V 2 2 /2Jc p -[V 2 2 - 2V 2 (V 2 - V)]/2Jc. 


-Oo* < 


T 2 + V 2 2/2Jc t 


£ dn 


1 + V 2 2/2Jc T 2 


V 


V 2 2 / 2 Jc 


T 2 + V 2 / 2Jc t 


V2® 


'2oo 


Jc To„, 

p 2oo 


-) 




= Jh£ 


t Jc T„ 

P 2® 


A - T 2o° A \t ' 

r u 5/ 


or, by equations (4.2) and (4.3), 


0 2 " 


Vo 2 

(H, 


L Jc P T 2co 


l) 

T 2oo 


■] 


* * 


(5.19) 


The term 0 2 /t', which by equation (5.19) is proportional to 
5 2 **/t’, is a small quantity of the first order } with the assumed order 
of approximation, equation (5.15) may be transformed into 


P*V, 


1 

0 Jc 


2o» 

a» Y 2°o|jc p T 2 J 


(Ho - 1) 


l 2tio 

^2o7 


8 2** “ p k V k 9 k 


or. 


1 

c 


<>a» T a. S 2 - ' (H * ' 1)T 


2oo 


&**=nVT *0 

2 p k V k ± 2» k 


(5.20) 
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The system of equations (5.18) and (5.20) gives the required 
system of equations for determining the two unknown magnitudes §.?** 

and Eh, as functions of the parameters of the boundary layer, of 
the external flow near the trailing edge of the airfoil, and of the 
density, velocity, and temperature at infinity behind the cascade. 

For the solution of this system of equations, it is noted that 
the unknowns are readily separated if equations (5.18) and (5.20) are 
divided one by the other. This procedure yields 


/VjA & V g^/jC p : 

VW *2-* 





or, when the Mach number Mg^ is introduced at infinity behind the 
cascade. 




.. V 2oo _ V 2 


OP _ 


V 2e 


'N 


a 2c» * 1 > J 5p T 2«> 


( 


itt 


V^_\ 2 2 (k - 1)M 2oo 2 - (h 2 - 1) /^2oo\ 

.. / 1^, ~ V-- O G TT / 




f(k - 1)M 2o 2 + 1 


1+ #*k 


/ (5.21) 




j 


This transcendental equation in Hg may be solved by one of the 
approximate methods in any concrete case. According to equation (5.20), 


'* „ 1 * - 1 ) M 2» 2 Plr y ir 

" 1 + (i - DM2,2 . a, p^v 2tt 


(5.22) 


The terms 5g** and H 2 referring to section 2, the location 
of which is unknown, are thus eliminated and expressed in terms of 
the magnitudes 5 k ** and 0^ either measurable or computable by 
any method of boundary- layer theory and in terms of the velocity, 
density, and temperature at the outer limit of the boundary layer 
near the trailing edge of the air foil and at infinity behind the 
cascade. The terms p' and R* may then be obtained with little 
difficulty by equations (4.8) and (4.9). 
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6. Approximate formulas for computing losses and resistances in 
cascade. - At these relatively small subsonic Mach numbers considered 
herein, the non! so thermal character of the flow in the wake behind 
the airfoil of the cascade can occur mainly through the heating or 
the cooling of the surface of the airfoil and not through the internal 
transformations of kinetic energy into heat. 

In order to verify this fact, equations (5.2l) is employed. 

The following notation is introduced for briefness: 




m 





> 


5 + H k 
2 




J 


( 6 . 1 ) 


Equation (5.2l), which is transcendental relative to a , then 
assumes the form 


m - £ _ 1 0 k / v 2o°^ 

1 + m 2 5 ^ \v k J 


(i+€ 


( 6 . 2 ) 


The unknown magnitude c is now expanded into a series in powers 
of the small parameter m. (For air the value of m at M 2qo < 0.7 
does not exceed O.l) 


e 


-0 


Cj_m 


+ £om 


(6.3) 


Substituting this series in equation (6.2) gives 


£ £ 0 + C 1 “ € l) m 



or 

“ € 0 + C 1 + € 0 " e l) m 
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By equating coefficients, the following equation is obtained: 



(6.5) 


for determining e^. Because of the assumption previously made on 
the small heat transfer from the surface of the airfoil in the cascade, 
the quantity is considered small for Mg^ = 0. The following 

equation, accurate to small quantities of the second order, is then 
obtained: 


“ *0 " 



E 


= A, 1 + € n In 


Vk 




( 6 . 6 ) 


where 



From equation (6.6), 


0 1 + A k In (V 2a /V k ) 


(6.7) 


( 6 . 8 ) 


The ratio generally differs little from unity* hence, 

the natural logarithm of this ratio is small so that Cq = - A k 
may be written without great error. 

By equating the coefficients of m to the first power, in 
equation (6.4), the expression e^=l + CQ=l-A k is obtained 

with the same degree of accuracy. 

The approximate equation is then 


e = - A k + (1 + A k )m « m - A k 


(6.9) 
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Equation (5.22) is now employed and in the new notation has 
the form 


s 2 ** - -2-L£- -2 & - 


( 6 . 10 ) 


According to equations (6.1l) and (6.7), 


r * * 

b 2 


i(5+Hj 




2oo 


V2J 


( 6 . 11 ) 


If it is assumed that at the trailing edge H k = 1.4, equations 
(6.11) and (6.7) assume the form 



( 6 . 12 ) 


From equation (4.8), an approximate formula for the losses is readily 
obtained: 



and therefore a corresponding approximate formula for the resistance 
differing from the right side of the previous equation only in the 
factor t. 
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The further possible simplifications of equation (6.13) are 
connected with the choice of devices for computing the characteristics 
of the boundary layer at the surface of the airfoil in the cascade 
and for taking into account the effect of the compressibility on the 
external flow. 


Translation by S. Reiss 
Rational Advisory Committee 
for Aeronautics 
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Figure 3. 
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